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We investigate impact of 7r°-production data at Relativistic Heavy Ion Collider (RHIC) and fu- 
ture E07-011 experiment for the structure function gi of the deuteron at the Thomas Jefferson 
National Accelerator Facility (JLab) on studies of nucleonic spin structure, especially on the polar- 
ized gluon distribution function. By global analyses of polarized lepton-nucleon scattering and the 
7i-°-production data, polarized parton distribution functions are determined and their uncertainties 
are estimated by the Hessian method. Two types of the gluon distribution function are investi- 
gated. One is a positive distribution and the other is a node-type distribution which changes sign 
at a: ~ 0.1. Although the RHIC tt" data seem to favor the node type for /S.g{x), it is difficult to 
determine a precise functional form from the current data. However, it is interesting to find that the 
gluon distribution l\g{x) is positive at large x (> 0.2) due to constraints from the scaling violation 
in gi and RHIC tt^ data. The JLab-E07-011 measurements for gf should be also able to reduce 
the gluon uncertainty, and the reduction is comparable to the one by RUN-5 vr^-production data 
at RHIC. The reduction is caused mainly by the error correlation between polarized antiquark and 
gluon distributions and by a next-to-leading-order (NLO) gluonic effect in the structure function 
gf. We find that the JLab-E07-011 data are accurate enough to probe the NLO gluonic term in gi. 
Both RHIC and JLab data contribute to better determination of the polarized gluon distribution 
in addition to improvement on polarized quark and antiquark distributions. 

PACS numbers: 13.60.Hb,13.88.-|-e 



I. INTRODUCTION 

Quark and gluon contributions to the nucleon spin are 
described by polarized parton distribution functions (po- 
larized PDFs) and their first moments. It became clear 
that only a small fraction of nucleon spin is carried by 
quarks and antiquarks. Therefore, a large gluon polar- 
ization or effects of orbital angular momenta should be 
possible sources for explaining the origin of the nucleon 
spin. The gluon polarization is expected to be clarified in 
the near future, whereas it would take time to determine 
the effects of the orbital angular momenta. 

Polarized PDFs have been investigated by global anal- 
yses of data on polarized lepton-nucleon deep inelastic 
scattering (DIS) and proton-proton collisions P, [3, H, 
0) S @|- Polarized quark distributions are determined 
relatively well; however, the polarized gluon distribu- 
tion A(7(x) is not accurately determined. Here, x is the 
Bjorken scaling variable, and Ag(a;) is the difference be- 
tween the gluon distribution with helicity parallel to that 
of parent nucleon and the one with helicity anti-parallel. 

The gluon distribution contributes to the structure 
function gi as a higher-order effect in the expansion by 
the running coupling constant Ug of quantum chromody- 
namics. The unpolarized gluon distribution has been de- 
termined primarily by the dependence of F2 at small 



X, where is defined by the momentum transfer q by 
= —q^ in lepton scattering. The kinematical range of 
X and is still limited for in determining ^g{x) by 
the scaling violation, so that the determination of ^g{x) 
is difficult from the scaling violation. Nonetheless, it 
is noteworthy that there are differences between the 
COMPASSS and HERMES data for gi in the range of 
X ^ 0.05. Such differences could be used for con- 
straininga gluon polarization at large x as pointed out 
in Ref. Q. However, this idea should be tested by fu- 
ture measurements on the scaling violation in gi because 
the differences could originate also from higher-twist 
effects 3. 

Other types of measurements are needed to improve 
the situation of /S.g{x). There were measurements on 
l^g{x) in lepton-nucleon scattering by observing high- 
Pt hadrons [3, H, and open-charm events fl^ . These 
data provided constraints on the gluon polarization at 
X ~ 0.1. They indicated that the ratio Ag{x)/g{x) is 
small although experimental errors are still large. Mea- 
surements at Relativistic Heavy Ion Collider (RHIC) are 
also important for constraining the gluon polarization. 
For example, 7r° and jet-production data [ll|, [l2| in po- 
larized proton-proton collisions are valuable for the de- 
termination of Ag{x). In fact, we showed that the tt^ 
data play an important role in reducing the uncertainty 



of Ag{x) by a global analysis including the vr*^ data in 
addition to the DIS data . Certain fragmentation func- 
tions are used in the analysis of 7r°-production processes. 
We should be careful that gluon and light-quark frag- 
mentation functions have large uncertainties at small 
or small pr ■ 

In order to determine the polarized PDFs including 
the gluon polarization, precise measurements are needed 
also in DIS. After our previous analysis new DIS 
data are reported by the CLAS 14], HERMES [H, and 
COMPASS [iBl collaborations. In future, the structure 
function gi for the deuteron will be accurately measured 
at the Thomas Jefferson National Accelerator Facility 
(JLab) by the proposed experiment E07-011 ■ The 

measurements at JLab should be valuable for reducing 
large uncertainties of the polarized antiquark and gluon 
distributions because they cover a large-x region at small 

In this work, we determine the polarized PDFs by 
global analyses of the data for spin asymmetries Ai in po- 
larized lepton-nucleon DIS and double spin asymmetries 
Ai^i^ of the TT*^ production in proton-proton collisions. In 
particular, we focus our discussions on determination of 
the polarized gluon distribution function. The purposes 
of this article are the following. 

(1) The RHIC 7r°-production data could indicate a 
node- type distribution for Ag(x), so that two types 
of distributions are studied in the global analysis by 
taking a positive distribution for Ag(x) and a node- 
type distribution which changes sign at x 0.1. 

(2) Among various future experimental projects, we in- 
vestigate impact on the polarized PDFs, especially 
on the polarized gluon distribution function, from 
precise measurements on the structure function gi 
of the deuteron by the JLab-E07-011 experiment 
[13, m. The E07-011 data will be so precise that 
Ag(x) could be constrained. 

(3) The polarized PDFs obtained with the E07-011 
data are compared with the ones including the vr"- 
production data especially in reducing uncertain- 
ties of Ag{x) and its first moment. 

We provide a new AAC08 (Asymmetry Analysis Collab- 
oration in 2008) library by the current global analyses 
because updated AAC distributions have not been pro- 
vided since the 2003 version [1, [l^ . 

This article is organized as follows. In Sec. |TT1 three 
data sets are introduced for investigating roles of the 
RHIC 7r° and JLab data, and our analysis method is 
explained. Analysis results are discussed in Sec. Illli and 
they are summarized in Sec. IIVI 



II. ANALYSIS METHOD 

In order to investigate roles of the RHIC tt*^ and E07- 
011 data, three data sets in Table U are prepared. The 



TABLE I: Used data sets in our global analyses. The "DIS" 
indicates current DIS data on gi . The RHIC tt" data are taken 
from Ref. The JLab-E07-011 data are calculated by 

using estimated errors in Refs. [TtI, [isj and set- A results, and 
actual vales are shown in Table HIl The notation Q indicates 
an included data set. 



Analysis set 


DIS 


RHIC tt" 


E07-011 


A 


O 






B 


o 


o 




C 


o 




o 



notation "DIS" indicates all the present DIS data on gi, 
and they are almost the same data used in the AAC06 
analysis [H]. There are minor changes. The same data 
are used for SLAC, EMC, SMC, and HERMES; however, 
updated ones are used for the deuteron by the COMPASS 
[l6| . In addition, the CLAS data are included in 
the DIS set. This is called set A; the one with 7r° data 
is set B; the one with the E07-011 is set C. The data 
for the RHIC 7r° and expected E07-011 are taken from 
Refs. and [13, [11, respectively. The RUN-5 data 
are used for the tt" production. Preliminary RUN-6 data 
by the PHENIX collaboration [2^ and jet data of the 
STAR collaboration are not included in our analyses. 
We may include these data in our future global analysis. 
Comparing three analysis results, we try to clarify roles 
of the RHIC 7r° and E07-011 data in the determination of 
the polarized PDFs, especially on Ag(x), in comparison 
with the results by the set A. 

The expected data for Af of the E07-011 experiment 
are estimated in the following way. The asymmetries 
are assumed to be the same as the ones of the analysis 
A: Af(a:;,(32)Eo7-oii = Af{x,Q'^)sct-A- Expected errors of 
the E07-011 experiment are estimated for {Sgf/ gf)EQ7-oii 
in the proposal [13, [3 , and they are converted to the 
errors of Af by 



dAf{x,Q^)E07-oii 



6g1ix,Q^ 



gfix,Q^) 



Ai{x,Q\,t- 



E07-011 



(1) 

Estimated values and their errors are shown for the asym- 
metry Af{x,Q^) in Table [Til It should be noted that 
there are two solutions for the polarized PDFs, the pos- 
itive and node-type distributions, so that there are two 
sets of Af data in Table HH 

The kinematical region of the used DIS data for pro- 
ton, neutron, and deuteron is shown in Fig. [T] by a; and 
Q2. The CERN data (EMC, SMC, COMPASS) cover a 
wide X region with relatively large Q^. The are especially 
important for determining the polarized PDFs at small x 
{x - 0.01). The SLAC (E130, E142, E143, E154, E155), 
DESY (HERMES), and JLab (Hall-A, CLAS, E07-011) 
data are at relatively large x with small Q^. The are es- 
pecially important for determining the polarized PDFs at 
large x and also possibly the polarized gluon distribution 
in comparison with the CERN data. 
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TABLE II: Expected asymmetries Af{x,Q^) and their errors 
in the E07-011 experiment. They are calculated by using 
the set- A analysis results and estimated errors Sgf/gf. The 
details are explained in the text. 



X 






Af 


6Af 


At 


SAf 




(GeV^) 




positive 


node 


0.175 


1.4 


0.020 


0.0841 


0.0017 


0.0819 


0.0016 


0.25 


1.9 


0.014 


0.1568 


0.0022 


0.1475 


0.0021 


0.35 


2.5 


0.010 


0.2522 


0.0025 


0.2453 


0.0025 


0.45 


3.0 


0.010 


0.3416 


0.0034 


0.3414 


0.0034 


0.55 


3.7 


0.015 


0.4308 


0.0065 


0.4384 


0.0066 



The details of our analysis method are explained in 
our previous article so that only a brief outline is 
explained in the following. The polarized PDFs are given 
at the initial scale (= Qq), where it is taken Qq = 1 
GeV^, with a number of parameters: 



Afix, Ql) = [Sx'' - ti{x- - x'^)]f{x, Ql) 



(2) 



where i5, k, v, and ^ are parameters to be determined 
by a analysis, and f{x) is the corresponding unpolar- 
ized PDF [2l| . Flavor symmetric antiquark distributions 
are assumed for the polarized PDFs at the initial scale 
although the unpolarized antiquark distributions are not 
flavor symmetric [22| . They are evolved to experimental 
points by the DGLAP (Dokshitzer-Gribov-Lipatov- 
AltarcUi-Parisi) equations [23j . 

The parameters are then determined by a analy- 
sis of experimental data on spin asymmetries so as to 
minimize 



E 



(3) 
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10: 



EMC-p 

SMC-p 

EI3D-P 

EI43-P 

ElSS-p 

HERMES-p 

CLAS-p 

E143-d 

ElSS-d 



H SMC-d 

HERMES-d 

X COMPASS-d 

□ CLAS-d 

X E07-011-d 

h HERMES-n 

=!■ E142-n 

+ ElS4-n 

X Hall-A-n 
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Here, the spin asymmetries are Ai in the DIS and AJ^j^ 
in the tt^ production. The analyses are done in the next- 
to-leading order (NLO) of the running coupling constant 
as, and the modified minimal subtraction (MS) scheme 
is used. The fragmentation functions of the HKNS07 
(Hirai, Kumano, Nagai, Sudoh in 2007) are used 
for describing the fragmentation into the pion. However, 
uncertainties of the fragmentation functions arc not in- 
cluded in this analysis. The uncertainties of the polarized 
PDFs are estimated by the Hessian method: 



[5F{x)r = Ax'Y. 



fdF{x,a) 



dF{x,a) 



0.001 



0.01 



0.1 



(4) 

where F{x) is a polarized PDF, Hij is the Hessian, the 
parameters are denoted (i=l, 2, N), and a indicates 
the parameter set at the minimum point. The value 
of Ax^ is taken as Ax^ = 12.65 so that the uncertainty 
indicates the one-cr-error range for the normal distribu- 
tion in the eleven-parameter space 0, [1] . The details of 
the analysis conditions and the uncertainty-estimation 
method are found in Ref. [51- 

We comment on our choice of Ax^- As explained in 
Ref. [l^l, Ax^ ~ N gives one-a error for a multivari- 
ate normal distribution of fit parameters. We took this 
Ax^ value for estimating the uncertainties by consid- 
ering it as an appropriate quantity for showing errors 
of overall functional behavior. However, this choice is 
not unique in global analyses of PDFs [l^, [2^. For 
example, Ax'^ = 1 is chosen in the COMPASS analy- 
sis for Ag{x) The difference between the choices, 
Ax^ ^ N and Ax^ = 1, comes from the fact that the 
likelihood function in the multiparameter space is dif- 
ferent from the likelihood function for a physics quan- 
tity, for example, Ai{x,Q^) at a fixed kinematical point 
of X and Q^. If the latter likelihood is plotted as a 
function of this physical quantity (one degree of free- 
dom) Ai{x,Q^), the one-fj (68%) range as the statistical 
error of the function corresponds to the Ax^ = 1 re- 
gion. However, our experience indicates that the choice 
Ax^ = 1 could be an underestimation of the uncertainties 
in comparison with experimental errors and their varia- 
tions (see Figs. 1 and 3 in Ref. d). Therefore, CTEQ 
(Coordinated Theoretical/Experimental Project on QCD 
Phenomenology and Tests of the Standard Model) and 
MRST (Martin-Roberts-Stirling-Thorne) use larger Ax^ 
values ( Ax^ = 50 or 100) [1^ as an effective way to show 
the error range. Our choice A^^ iV is a conservative 
estimation of the uncertainty range by considering error 
correlations among the parameters. If one would like to 
use the choice A^^ = 1, one may scale down our polar- 
ized PDF uncertainties by l/y/Ax^. 



FIG. 1: Kinematical region of the DIS data is shown by x and 
Q^. The notations p, n, and d indicate proton, neutron, and 
deuteron, respectively. 



III. RESULTS 

Three types of analyses are done with the different 
data sets in Table HI The x^ values are shown for each 
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TABLE III: Numbers of data and values are listed for the three analysis sets with the data in Table U] The positive and 
node indicate a positive distribution and a node-type one, respectively, for Ag{x). The notations p, n, and d indicate proton, 
neutron, and deuteron, respectively. 



Data set 


No. of 
















data 


Set A 




Set B 




Set C 








Positive 


Node 


Positive 


Node 


Positive 


Node 


EMC (p) 


10 


5.17 


4.52 


5.18 


4.57 


5.17 


4.52 


SMC (p) 


59 


56.34 


53.79 


56.77 


53.95 


56.34 


53.79 


E130 (p) 


8 


5.13 


4.84 


5.23 


4.97 


5.13 


4.84 


E143 (p) 


81 


59.97 


59.78 


60.43 


60.30 


59.97 


59.78 


E155 (p) 


24 


33.07 


27.38 


31.76 


26.19 


33.06 


27.39 


HERMES (p) 


9 


3.58 


3.69 


3.59 


3.68 


3.58 


3.69 


CLAS (p) 


10 


14.34 


15.11 


17.81 


21.65 


14.32 


15.10 


SMC (d) 


65 


57.39 


56.65 


57.82 


56.86 


57.39 


56.65 


E143 (d) 


81 


91.15 


91.86 


89.84 


89.85 


91.16 


91.87 


E155 (d) 


24 


18.60 


22.44 


18.53 


23.32 


18.60 


22.44 


HERMES (d) 


9 


11.19 


8.10 


11.93 


7.65 


11.20 


8.10 


COMPASS (d) 


15 


8.73 


10.49 


9.30 


11.89 


8.73 


10.49 


CLAS (d) 


10 


5.73 


4.76 


5.55 


5.53 


5.73 


4.76 


E07-011 (d) 


5 










0.00 


0.00 


E142 (n) 


8 


2.60 


2.48 


2.78 


2.68 


2.60 


2.48 


E154 (n) 


11 


3.32 


4.09 


3.07 


3.12 


3.32 


4.10 


HERMES (n) 


9 


2.26 


2.41 


2.18 


2.36 


2.26 


2.41 


Hall-A (n) 


3 


3.09 


3.49 


2.89 


2.70 


3.09 


3.48 


DIS total 


441 


381.66 


375.87 


384.65 


381.27 


381.66 


375.87 


PHENIX (tt") 


10 






12.43 


11.32 






Total 


451 


381.66 


375.87 


397.08 


392.60 


381.66 


375.87 


(xVd.o.f.) 




(0.90) 


(0.88) 


(0.91) 


(0.90) 


(0.89) 


(0.87) 



TABLE IV: First moment of the polarized gluon distribution function and quark spin content AS — jj^ dx [Aqi{x) + Aqi{x)] 
at Q2 = 1 GeV^ 





Set A 


Set B 


Set C 




Positive Node 


Positive Node 


Positive Node 


AE 


0.24 ± 0.07 0.22 ± 0.08 


0.26 ± 0.06 0.25 ± 0.07 


0.24 ± 0.05 0.22 ± 0.05 


AG 


0.63 ± 0.81 0.94 ± 1.66 


0.40 ± 0.28 -0.12 ± 1.78 


0.63 ± 0.45 0.94 ± 1.09 



analysis in Table [TTTl The total values per degrees of 
freedom are in the range x^/d.o.f.=0.87— 0.91, so that all 
the analyses are successful in explaining the data. The 
value for the expected E07-011 data in the analysis 
C is very small (~ 10"'') because the values of the data 



TABLE V: The first moment AG and its uncertainty SAG 
at = 1 GeV^ in the range 0.1 < a; < 1. Here, the inte- 
grals are calculated in the limited x range by AG(a; > 0.1) = 
fl^dxAg{x). 



Function 


Set 


AG(x > 0.1) 


5AG{x > 0.1) 


sAa(x>o.i) 

AG(x>0.1) 


Positive 


A 


0.53 


0.72 


1.36 




B 


0.36 


0.26 


0.71 




C 


0.53 


0.38 


0.73 


Node 


A 


0.87 


0.89 


1.02 




B 


0.40 


0.31 


0.77 




C 


0.87 


0.47 


0.54 



are taken from the analysis-A results and there are no 
statistical variations. The polarized PDFs of the anal- 
ysis C should be close to the ones of the analysis A. 
Comparison between A and C results is useful for dis- 
cussing constraints from the E07-011 data on ^g{x) de- 
termination by noting differences between their uncer- 
tainties. First moments of the polarized PDFs are listed 
in Table ITVl The first moments of Auv{x) and Ady{x) 
are fixed by semileptonic decays: dx Auy{x) = 0.926 
and Jq dx Ady{x) ~ —0.341 4]. In order to discuss con- 
straints from the RHIC-tt" and JLab-E07-011 data on 
the gluon distribution at relatively large x, we also show 
the "first moment" calculated by the integral AG{x > 
0.1) = J^^dxg{x) in TableEl 



A. Impact of RHIC tt" data 

First, effects of the RHIC tt*' data [ll| are shown in 
comparison with the polarized PDFs obtained only by 
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the DIS data. Using the prehminary PHENIX 7r° data 
of the RHIC RUN-5, we have already shown their impact 
on the gluon distribution in Ref. [5|. The current anal- 
ysis uses the published RUN-5 data. Analysis-B results 
are compared with analysis-A ones in Fig. [21 where the 
polarized PDFs xAuy, xAdy, xAq, and xAg and their 
uncertainties are shown at Q^ = l GeV^. The set A uses 
only the DIS data and the set B also includes the tt^ data 
as listed in Table HIl 

The polarized valence-quark distributions Au„ (x) and 
Adv{x) are not changed even if the RHIC 7r° data are 
included in the analysis. However, there are significant 
effects on the other distributions, especially on the po- 
larized gluon distribution. The gluon distribution has a 
huge uncertainty if it is determined only by the DIS data 
(analysis A); however, the uncertainty band becomes sig- 
nificantly small in the analysis B because of the addition 
of the 7r° data. It is known that gluon-gluon interaction 
subprocesses dominate the Tr^-production cross section 
especially at small pr, so that the 7r° data are sensitive to 
the gluon distribution. We notice that the antiquark un- 
certainty is slightly reduced in the set B and it is caused 
by the error correlation between the antiquark and gluon 
distributions as pointed out in Ref. 

The double spin asymmetry for the tt" production in- 
dicates a negative value at pr = 2.38 GeV. It suggests 
a node type distribution which vanishes at x ~ 0.1 as 
investigated in Refs. 0, d, HI]. We also made a global 
analysis by assigning an initial distribution to a node- 
type one for Ag{x). Minimum values per degrees are 
slightly smaller in the node-type results than the values 
in the positive ones. In the node-type results of the anal- 
ysis A, the gluon distribution is positive at a; > 0.2 and 
is slightly negative at small x. The negative gluon dis- 
tribution at X < 0.1 becomes smaller in the node- type of 
the analysis B. 




-0.5 
0.01 



S -O-Ol 
-0.02 
-0.03 



§- Set A (DIS) 
Set B (DIS + K*) 



Ag( positive) 




As noticed in Ref. [5], the differences between HER- 
MES (Q2 ^ 1 GeV2) and COMPASS (Q^ ~ 6 GeV^) 
data of g'l at X r-^ 0.05 could suggest a positive Ag{x) 
at large x (> 0.2) although it may be explained by a 
higher- twist effect Q . The gluon distribution contributes 
to gi (x) as a NLO effect and the differences could be ex- 
plained by this NLO effect with Ag{x) > at a; > 0.2. In 
fact, all our global analyses indicate positive distributions 
at large x. 

The node-type distributions are compared with the 
positive ones at = 1 GeV^ in Fig. [31 Both distri- 
butions are determined by the set-B analysis. It is clear 
from this figure that the gluon distribution at x < 0.1 
cannot be determined by the current data including the 
RHIC tt'^ data. In fact, the uncertainty of AG is large in 
the node-type distribution and it is obvious that a large 
contribution comes from the small x region. Because 
there is no constraint at small x, uncertainty reduction 
of AG is not very obvious due to the tt" data from the 
set A to set B as shown in Table IIVI 

In order to investigate a possible improvement on AG, 
we discuss a contribution to the first moment from the 
region of x > 0.1 in Table fVl The relative uncertainties 
become smaller in the set B from the set A: SAG{x > 
0.1)/AG{x > 0.1) = 1.36 0.71 in the positive dis- 
tribution and 1.02 0.77 in the node one. There are 
significant reductions in the uncertainties (24%~47%), 
which suggests that the 7r° data impose a significant con- 
straint in the gluon distribution at x > 0.1 in addition 
to the aforementioned DIS constraint from the scaling 
violation. 

We discuss a comparison with our previous work of 



Set B (Positive) 
Set B (Node) 



< 



1 






= 1 GeV^ 





FIG. 2: Comparison of the polarized PDFs between the anal- 
yses A and B. They are shown by the dashed and solid curves 
for the analyses A and B, respectively, at — 1 GeV'^. Their 
uncertainties are shown by the shaded bands. 



FIG. 3: Polarized antiquark and gluon distributions by the 
analysis B which includes the RHIC yr^-production data. The 
distributions and their uncertainties are calculated at = 1 
GeV^. Positive and node- type solutions are shown by the 
solid and dashed-dot curves, respectively. 
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AAC06 ^ . The positive distributions of the set B and 
their uncertainties are compared with the AAC06 results 
in Fig. m There are differences between the used data 
sets. As explained in Sec. HI additional data, which are 
not used in the AAC06 analysis, are those of the CLAS 
[i3 and COMPASS [3 collaborations. The PHENIX 
7r° data [HI were also slightly changed from the prelim- 
inary ones at the stage of the AAC06 publication. The 
used fragmentation functions are also changed from the 
functions by KKP (Kniehl, Kramer, and Potter) for the 
ones by the HKNS07. 

As shown in Fig. [H the distributions and their un- 
certainties are almost the same. There are slight im- 
provements in the valence-quark and antiquark distribu- 
tions, which is manly due to the accurate measurements 
at medium x by the CLAS and COMPASS collabora- 
tions. The slight reduction in the Ag{x) uncertainty 
band should come from the error correlation between the 
quark and gluon distributions. In addition, there are ef- 
fects due to changes in the tt" data from the preliminary 
version and in the fragmentation functions. However, 
these effects are not significant because the distributions 
themselves are almost the same. We also compared the 
AAC06 and current distributions for the node type, but 
they are also almost the same in both distributions and 
uncertainties. 

It should be noted that the description of a 7r° pro- 
duction cross section depends much on pion fragmen- 
tation functions. Here, the HKNS07 functions [1^ are 
employed in our global analyses for describing the pion 
cross section. As noted in Ref. ^1^, there are large un- 
certainties in the fragmentation functions, especially in 
the gluon and light-quark functions. Effects of such un- 
certainties are shown in Fig. [5] together with the RHIC 
Tr'^-production data and the asymmetries calculated by 
using the fragmentation functions of Kretzer and KKP 



- AAC06 

- Set B (DIS + 



Ag( positive) 



0.06 






Q = 1 GeV' 



FIG. 4: Comparison of the polarized PDFs between the cur- 
rent analysis B and AAC06 for the positive Ag{x). They are 
shown by the dashed and solid curves for the AAC06 and the 
current analysis (AAC08), respectively, at — 1 GeV^. Their 
uncertainties are shown by the shaded bands. 
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FIG. 5: Set-B results are compared with the PHENIX data 
in the longitudinal double spin asymmetry oi p + p ^ tt'^ + X 
[nj . The solid curves and the shaded bands indicate theo- 
retical asymmetries and their uncertainties by the HKNS07 
fragmentation functions. The results by the KKP and Kretzer 
fragmentation functions are shown by the dashed and dashed- 
dot curves, respectively. 



as other examples. 

The pion production is sensitive to the gluon fragmen- 
tation function, which has a large uncertainty because it 
has been determined mainly by the scaling violation of 
e^"e~ h + X data. However, we expect that the situ- 
ation will improve if the Belle and Babar collaborations 
provide low-energy e"'"e~ data for clarifying the scaling 
violation. In our analysis of the pion production, the 
NLO effects are approximately taken into account as a 
if-factor by first calculating the leading-order (LO) cross 
sections. Therefore, the LO fragmentation functions are 
used in the analysis, and the asymmetry has the large 
uncertainty bands as shown in the figure. The gluon frag- 
mentation function is determined more precisely in the 
NLO [l^, so that the uncertainties of AJ^j^ are reduced 
about 1/3. In any case, there are significant effects on the 
spin asymmetry AJ^j^ from the fragmentation functions, 
which indicates the importance to access the uncertain- 
ties from the fragmentation functions in determining the 
error of Ag(x). 



B. Impact of JLab-E07-011 data 

We discuss the impact of the proposed experiment 
JLab-E07-011 [l3| to measure the structure function gi 
for the deuteron on the determination of the polarized 
PDFs. The expected E07-011 data [H in Table HI are 
included in the analysis C. In order to find an improve- 
ment to the analysis A with the current DIS data, the 
polarized PDFs of the analyses C are compared with the 
analysis-A results in Fig. [6] at Q^=l GeV^. The uncer- 
tainties are shown by the shaded bands. 

Both distributions are almost the same because the 
E07-011 data are assumed to be the same asymmetries 
obtained in the analysis A. The uncertainties are sim- 
ilar in the polarized valence-quark distributions in Fig. 
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[SI whereas the uncertainty of the antiquark distribution 
becomes smaher in the analysis C than the one of the 
analysis A. It is also clear that the uncertainty of the 
gluon distribution is significantly reduced in the analysis 
C. 

There is a possibility that the reduction of the gluon 
uncertainty due to the E07-011 data could come from the 
dependence in comparison with other data. In order 
to discuss such a possibility, we show the dependence 
of within the x range of the E07-011 data in Fig. [71 
Two X ranges are shown. One is at 0.166 < a; < 0.182 
and the other is at 0.416 < x < 0.50. The solid and 
dashed curves are theoretical set-C results for the positive 
Ag(a;) distribution and an artificial one with ^gix) — 
at the initial scale Q\. We notice that differences be- 
tween these curves are very small in comparison with the 
experimental errors. The results suggest that the current 
data should not be accurate enough to probe the polar- 
ized gluon distribution from the dependence in the x 
region, 0.1 < a; < 0.5. Therefore, the uncertainty reduc- 
tion in due to the E07-011 data should come from 
other sources. 

There are two possibilities for the reduction of the 
gluon uncertainty. First, it could be due to the error 
correlation between the antiquark and gluon distribu- 
tions. A more accurate determination of antiquark dis- 
tributions results in the improvement on the gluon deter- 
mination through the error correlation between polarized 
antiquark and gluon distributions. 

Second, it could be due to the NLO term with the 
gluon distribution. In Fig. [51 wc show the ratio of the 
gluon NLO term to g\: 



1 1 



5i(x,Q2)2^ 
where is the quark charge of the flavor i, and ACg is a 



gluonic coefficient function. The ratios are calculated by 
using the polarized PDFs of the analysis C, and they are 
shown by solid boxes and triangles. The NLO corrections 
are calculated at the corresponding experimental x and 
values. The experimental errors are shown by the 
ratio 



51 (a;, 



(6) 



where Q^) and Q^) are an experimental 

spin asymmetry and its error. 

It is obvious from Fig. [51 that the gluonic NLO ef- 
fects are within the experimental errors of the current 
CLAS data lj| for both proton and deuteron. How- 
ever, the proposed experiment E07-011 is very accurate, 
and expected errors are much smaller than the NLO ef- 
fects shown by the boxes and triangles. It suggests that 
the E07-011 measurements should be valuable for deter- 
mining not only the quark and antiquark distributions 
but also the gluon distribution. However, higher-twist 
effects become apparent at small , so that careful con- 
sideration is needed also for them as well as the gluonic 
NLO contributions if such high-precision data are taken 
at JLab. 

Another global analysis is made by removing only the 
CLAS data from the data set A, and results indicate that 
the I^g{x) uncertainty is slightly reduced. From Fig. [51 
we find that such a reduction due to the CLAS data 
is caused mainly by the error correlation between anti- 
quark and gluon distributions because the NLO effects 
are within the experimental errors. On the other hand, 
the reduction due to the E07-011 data in Fig. [6lis caused 
by both the error correlation and the NLO contributions 
because the NLO effects are much larger than the errors. 
It is important to find that the E07-011 data are very 
accurate to probe the higher-order gluonic effects in the 



Q^= 1 GeV^ 




Ag(positive) 





SetA(DlS) 
— SetC(DIS + E07-011) 



Ag( positive) 




FIG. 6: Comparison of the polarized PDFs between the anal- 
yses A and C. They are shown by the dashed and solid curves 
for the analyses A and C, respectively, at Q^ — \ GeV'^. Their 
uncertainties are shown by the shaded bands. 
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FIG. 7: dependence is shown in the spin asymmetry Ai 
for the deuteron. The data at 0.166 < a; < 0.182 and 0.416 < 
X < 0.50 are shown in comparison with the theoretical curves 
of the analysis C. The solid curves are calculated by the set- 
C distributions with Ag{x) > 0, and the dashed ones are 
by the set-C distributions with a different gluon distribution, 
Ag{x) = at Qq — 1 GeV^. The theoretical curves are 
calculated at a; = 0.175 and x — 0.45. 
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structure function gi. 

We comment whether or not such a gluonic NLO term 
might be effectively absorbed into variations of polarized 
antiquark distributions because it is generally considered 
to be difhcult to separate both distributions by methods 
other than the scaling violation. However, it is not the 
case in our analysis with the E07-011 data. As clearly 
shown in Fig. [6l the reduction of the gluon uncertainty is 
much larger than the antiquark one. It cannot be simply 
explained by the error correlation between the antiquark 
and gluon distributions. It suggests that the E07-011 
data should have a significant impact on the gluon dis- 
tribution and that the effect should not be absorbed into 
the antiquark distributions. 

Next, we compare the gluon distributions A.g{x) of the 
analysis C with the analysis-B distributions in Fig. [S] 
This figure is intended to show the impact of the E07- 
011 data in comparison with effects of the RHIC 7r° data. 
We found in Fig. [9] that both B and C distributions are 
similar. However, it should be noted that the positivity 
condition |A5(a;)| < g{x) is not satisfied for the node- 
type solution in the analyses B and C at large x. It was 
difficult to get a converging result within the positivity 
condition especially for the parameter 6g, so that the 
condition is not imposed in the analysis. We also notice 
that the relative uncertainties for C are roughly the same 
with the ones for B. This fact indicates that the effects 
of the E07-011 data on the determination of ^g{x) are 
comparable with those of the 7r° data. 

Mechanisms for the reductions of the gluon uncertain- 
ties are completely different in both cases. In the analysis 



B, the reduction is due to the gluon-gluon scattering sub- 
processes in the tt*^ production, and it is due to the error 
correlation and the NLO gluon term of gi in the analy- 
sis C. Nonetheless, it is noteworthy that high-precision 
lepton-scattering data can be used for determining the 
polarized gluon distribution and that the improvement 
is as good as the current RHIC 7r° data. 

We comment on possible higher-twist effects. Because 
the JLab data are taken at relatively small Q^, higher- 
twist corrections may affect the results. Such effects are 
not included in our leading- twist analysis. There were 
studies on the higher-twist effects in the structure func- 
tion gi ^, j|] . We also have done a preliminary analysis on 
such effects; however, they are not precisely determined. 
Namely, uncertainties of the higher-twist corrections are 
large, which indicates that they may not be determined 
by the current experimental data as pointed out in Ref. 
[3] ■ If additional parameters are introduced in the higher- 
twist term, they have very large errors, which makes it 
difficult to get a reliable polarized PDF set by the current 
global analysis. We will work on an analysis for obtaining 
a possible higher-twist correction from the experimental 
data. 

Finally, we discuss the first moment of the polarized 
gluon distribution, AG (= Jp Ag{x)dx). The first mo- 
ments and their uncertainties are shown in Table IIVI for 
all the analyses. We find that the uncertainties become 
smaller in the set C from the set A: 5AG^ 0.81 — > 0.45 
in the positive distribution and 1.66 — > 1.09 in the node 
one. In order to investigate a possible constraint on AG 
in the large-x region, we also showed the first moment 
in the kincmatical region of x > 0.1 in Table fVl The 
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FIG. 8: The ratio of the gluon NLO-correction term to the 
polarized structure function gi{x,Q^) in Eq. ([5]). It is com- 
pared with experimental errors, which are shown by the ratio 
5giix,Q^)/giix,Q^) in Eq. ©. 



FIG. 9: Polarized gluon distributions and their uncertainties 
of the analysis C are compared with the ones of the analysis 
B at — 1 GeV^. Both positive and node-type distributions 
are shown. The set B and C results are shown by the solid 
and dashed curves. 
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TABLE VI: The first moment AG and its uncertainty SAG 
at = 1 GeV^ in the range 0.1 < a; < 1. The notations are 
the same as the ones in Table [V] 



Function 


Set AG{x > 0.1) 


SAG{x > 0.1) 


d/\U(x>0.1) 
AG(x>0.1) 


Positive 


C 0.49 


0.47 


0.96 


Node 


C" 0.74 


0.52 


0.70 



relative uncertainties become smaller in the set C from 
the set A: SAG{x > 0.1)/AG(x > 0.1) = 1.36 ^ 0.73 
(0.71 in set B) in the positive distribution and 1.02 
0.54 (0.77 in set B) in the node one. 

The analysis C has been made with the projected E07- 
011 data which are created so as to agree with the set-A 
spin asymmetries. In actual experimental measurements, 
statistical fluctuations exist. Therefore, the reduction of 
the gluon uncertainty could be overestimated. In order to 
investigate such a possibility, we tried the following anal- 
yses by generating the E07-011 data with two different 
assumptions: 

• Analysis C: Projected E07-011 data are obtained 
by using the node-type distributions in the set A, 
and then a global analysis is made with positive- 
type initial distributions. 

• Analysis C": Projected E07-011 data are obtained 
by using the positive-type distributions in the set 
A, and then a global analysis is made with node- 
type initial distributions. 

Analysis results of C indicate that the polarized 
PDFs and their uncertainties are almost the same as 
the positive-type distributions of the analysis C in the 
valence-quark and antiquark parts. Determined PDFs 
of C" are also almost the same as the node-type ones 
of the set C. However, the gluon distributions and their 
uncertainties are slightly modified although x-dependent 
functional forms are not changed. The gluon first mo- 
ments and their uncertainties are shown in Table IVll As 
expected, the uncertainties become larger than the ones 
for the analysis C in Table El SAG{x > 0.1)/AG(a; > 
0.1) = 0.96 [0.73 in set C; 0.71 (B); 1.36 (A)] for the 
positive type; 0.70 [0.54 in set C; 0.77 (B); 1.02 (A)] for 
the node type. According to the C and C" results, the 
improvement from the set A is not as large as the one in 
the set C. However, it is still comparable with the RHIC 
improvement in B. 

It is important to find that the large uncertainty of AG 
is significantly reduced by the proposed E07-011 exper- 
iment according to the analysis C (and C, C"). The 
improvement of AG is especially clear in the integral 
over the region x > 0.1. Therefore, the proposed JLab 
experiment is valuable for the determination of AG in 
addition to the quark and antiquark moments. If the po- 
larized lepton scattering data are accurate enough, they 
can contribute to a better determination of the gluon spin 



content. It is interesting to find that the analysis-C un- 
certainties are almost the same as the ones in the analysis 
B. However, it should be noted that the uncertainties of 
the pion fragmentation functions and scale uncertainties 
are not taken into account in the error estimate of the 
analysis B. 

From these comparisons, we found that the JLab mea- 
surements of the the proposed E07-011 can contribute 
to reducing the uncertainty of the gluon spin content by 
about 47% in the analysis C (30% in C and C"). The 
RHIC TT*^ measurements also reduce the uncertainty, but 
the inclusive lepton scattering data are valuable in the 
sense that they are free from the ambiguities of the frag- 
mentation functions [Tsj . 

We should, however, mention that future RHIC exper- 
iments will also improve the determination of Ag{x) by 
measurements of pion, jet, and direct-photon production 
processes. To be fair with the RHIC-Spin and other high- 
energy spin projects, it is desirable that a global analysis 
should be made together with other expected data by the 
time of the future JLab-E07-011 data. Furthermore, the 
current estimation would indicate a slightly better de- 
termination of Ag(x) from the JLab data because some 
fluctuations are expected in actual measurements from 
the cental spin asymmetries. However, we believe that 
our analyses with the expected errors are good enough in 
comparing the obtained relative errors for the gluon spin 
content from the DIS and collider data. In this article, 
we simply showed the role of accurate DIS data, by tak- 
ing the JLab-E07-011 data as an example, in comparison 
with the current status because all the expected data are 
not available for our analysis and also it is not easy to 
judge which data will be taken before the JLab-E07-011 
data. 

Because it is one of our purposes to provide a polarized 
PDF library, which has not been provided since 2003, we 
did not include preliminary data such as the RHIC run-6 
data in our global analyses. In addition, semi-inclusive 
data are also not included because we focused our anal- 
ysis on the determination of Ag{x) and the results de- 
pend much on used fragmentation functions. In future, 
we expect to make a more complete analysis including 
the RHIC run-6 and semi-inclusive data. 



IV. SUMMARY 

We have investigated the impacts of 7r°-production 
data in polarized pp collisions at RHIC and measure- 
ments of gf in the future JLab experiment E07-011 on the 
determination of the polarized PDFs, especially the po- 
larized gluon distribution. Global analyses of the polar- 
ized DIS were done for determining the polarized PDFs 
and their uncertainties by using the DIS data on gi, the 
7r° data of the RHIC RUN-5, and the expected JLab- 
E07-011 data on gf. 

The RHIC 7r° data indicated a possibility of a node- 
type distribution for Ag{x) which changes sign at x ~ 
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0.1. Accordingly, we made two types of analyses with 
positive and node-type distributions for Ag(x). Simi- 
lar Xmin values were obtained by both analyses although 
the node- type distributions have slightly smaller Xmin- It 
suggests that it is difficult to determine a precise gluon 
distribution from the current data. However, it is inter- 
esting to find in our analyses that the distribution A.g{x) 
is positive at x > 0.2 because of scaling violation in gi 
and Tr'^-production data. We found a large uncertainty 
for Ag{x) at small x (< 0.1). The uncertainty of the po- 
larized gluon polarization is significantly reduced by the 
RHIC 7r° data. In particular, the 7r° data constrain the 
gluon moment calculated in the large-x region (x > 0.1), 
and the reduction of the gluon uncertainty is 2A%^A7% 
according to our analysis. 

There is a clear improvement in the determination of 
A.g{x) by the future measurements of the E07-011 experi- 
ment. It is partially due to the error correlation between 
the polarized antiquark and gluon distributions. How- 
ever, the major improvement comes from the fact that 
the measurements are accurate enough to probe the po- 
larized gluon distribution in the structure function 171. 
Namely, the experimental errors are much smaller than 
the typical gluonic NLO term in gi. The E07-011 data 
can reduce the uncertainty of the first moment AG by 
about 30—47% if only the inclusive DIS data are used 
in the analysis. This reduction is comparable to the ef- 
fect of the RHIC RUN-5 tt" data. However, the JLab 



E07-011 data could have an advantage over the hadron- 
production data because there is no uncertainty com- 
ing from the fragmentation functions and the choice of 
the hard scale. A better determination of Ag{x) should 
be also made by future measurements of pion, jet, and 
direct-photon production processes at RHIC. 

Precise DIS and collider measurements are important 
for determining the polarized PDFs, especially in the po- 
larized gluon distribution. From the analyses B, we pro- 
vide a code for calculating two sets of optimum polarized 
PDFs. The set-1 is for the positive-type PDFs in the 
analysis B, and the set-2 is for the node-type ones also 
in the analysis B. The code is provided at our web home 
page [ii]. 
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